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Abstract 

in _ 

£^ . For a tight binding model with nearest neighbour attraction and a small 

. orthorhombic distortion, we find a phase diagram for the gap at zero tem- 

in 

ON . perature which includes three distinct regions as a function of filling. In the 



first, the gap is a mixture of mainly <i-wave with a smaller extended s-wave 



S. 

. part. This is followed by a region in which there is a rapid increase in the 

c : 

s-wave part accompanied by a rapid increase in relative phase between s and 



d from to ir. Finally, there is a region of dominant s with a mixture of d 
and zero phase. In the mixed region with a finite phase, the s-wave part of 
the gap can show a sudden increase with decreasing temperature accompa- 
nied with a rapid increase in phase which shows many of the characteristics 
measured in the angular resolved photoemission experiments of Ma et al. in 
Bi2Sr2CaCu208 
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Many experiments in the last few years have been interpreted in terms of a gap with 
d x 2_ y 2 symmetry in the Cu02 planes. Among these is the penetration depth experiment of 
Hardy et al.0 which show a linear dependence over a large temperature range down to the 
lowest temperatures measured. This is illustrative of a large class of experiments indicating 
that the gap goes to zero on the Fermi surface but which are not sensitive to its phase. 
Many new experiments have been devised specifically to test the phasei~i of the gap over 
the Brillouin zone with mixed results. Some favour d x 2_ y 2 symmetry, others extended s- 
wave. Perhaps the most detailed information on the momentum dependence of the gap 
has come from a variety of angular resolved photoemission data. Recently, Ding et al.H2l 
have found that the zero's in the gap are displaced somewhat from the main diagonals to 
either side of the (ir, tt) direction. This can arise if the gap is not pure d x 2_ y 2 but has an 
additional small s-wave component with zero relative phase between the two. This can 
be achieved if a small orthorhombic distortion is introduced in the Cu02 planeEUll and in 
addition the bilayer natureSS of B^S^CaC^Og is taken into account to give nodes on 
both sides of the diagonal. In the experiment of Ma et al.,0 the temperature variation of 
the gap is investigated in different directions in the Brillouin zone. In particular, the gap is 
found small and perhaps zero near T c in the (jr, 7r) direction and to increase rapidly below 
a reduced temperature of T/T c = 0.8. In the (0,7r) direction, the gap is large and follows 
more closely a BCS behaviour. 

In this paper, we investigate the phase diagram for simple BCS model with first nearest 
neighbour attraction with a view at understanding its phase diagram as a function of filling 
n and of making connections with the experimental situation described above. For the band 
structure, we employ a tight binding model with first and second nearest neighbour hopping. 
For tetragonal symmetry, we find at zero temperature, as a function of decreasing filling a 
region of pure d x 2_ y 2 symmetry which is followed by one of mixed d x 2_ y 2 plus an extended 
s x 2 +y 2 part with relative phase of 7r/2 between the two pieces. Finally at yet lower filling, 
there is a region of pure s x 2 +y 2 gap. If some small but finite orthorhombic distortion is 
introduced into the CuC>2 plane, the situation changes considerably in the sense that no 



pure d x 2_ y 2 region remains. Instead, it is replaced by a mixture of mainly d x 2_ y 2 but with a 
small admixture of s x 2 +y 2 with no relative phase so that in this region, the nodes in the gap 
are shifted somewhat off the (n, it) direction. As the filling is lowered at some particular 
critical value of n, the relative phase starts growing sharply from 0. In this regime there are 
only point nodes in the gap at (±7r/2, ±7r/2) in the Brillouin zone of the CuC>2 plane. This 
increase in relative phase is accompanied by a rapid increase in the s x 2 +y 2 part and similar 
sharp decrease in d x 2_ y 2. When the phase has increased to n, the s-wave part dominates 
and the gap will again have zeros on the Fermi surface. Adding an on site repulsion to 
the model changes the range of the various regions somewhat but more distinctively adds 
a small constant component to the gap which can be of the same or opposite sign as the 
s x 2 +y 2 part. 

The BCS gap equation determines the energy gap or order parameter Ak as a function of 
momentum k in the reciprocal lattice. To describe electrons in a 2-dimensional CuC>2 plane, 
we will use a tight binding dispersion relation with up to second nearest neighbours of the 
form 



with t the first nearest neighbour hopping parameter which is often assumed to be lOOmeV, 
B second neighbour hopping in units of t often taken to be B = 0.45 for YxE^CusOy 
and fj, is the chemical potential. The final parameter 5 in equation (1) is an orthorhombic 
distortion introduced into the CuC>2 plane which could be thought of as roughly modeling 
the existence of chains in YBaCuO. It has almost become conventional to eliminate the 
chemical potential in Eq. [I] in favour of the filling factor n defined as 



where the sum over k is over the first Brillouin zone of the CuC>2 plane and M is a normal- 
ization factor so that half filling corresponds to n — 1. 



£k = — 2i [cos(k x a) + (1 + 5) cos(k y a) 




(1) 




(2) 
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In momentum space, the BCS gap equation is given in terms of the pairing potential 
Vk_k' and the temperature [3~ l = ksT by 
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with E — Jsl + A k the quasiparticle energies in the superconducting state. 
The nearest neighbour part of the pairing potential has the form: 

= 2g x cos(k x - k' x ) + 2g y cos(k y - k' y ), (4) 

where a x and a y , the lattice parameters, are both taken to be a. In Eqs. |l] and [| it can be set 
equal to one if the momentum is measured in units of 27r/a. In Eq. |], g x and g y are coupling 
constant taken to be equal and their value set by the desired value of critical temperature 
T c or of gap amplitude. With the form in Eq. |] for the pairing potential, the BCS equation, 
Eq. H can be solved by convolution using fast Fourier transform techniques using 

A k = 77 E ^k-k'iv = T- 1 \T [V v ] ■ T [F v ]] , (5) 
■/v k , 

where J-\x\^ = Sk e lkr £k is the Fourier transform operator and = Ak tanh ((3Ek/2) jlE^. 

Another technique for solving the gap equation, Eq. |3], with the pairing potential of Eq. |] 
is to expand the interaction in terms of separable functions: 

= 2g x {cos(k x ) cos(£4) - sin(k x ) sin(^,)} 
+ 2g y {cos(k y ) cos(k' y ) - sm(k y ) sm(k' y )}. 

(6) 

The sin(k XtV ) terms do not contribute for singlet pairing. For a general separable inter- 
action of the form: 

vtt' } = E (7) 

n,m 

of which Eq. M is a special case, we can write Eq. as: 
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A k = E < E V n,mYf E VvFv, (8) 

and we see that Ak has no more components rfc than there are in the interaction. We will 
assume the rf£ to be an orthonormal set and only a few are needed to represent Eq. |[ The 
fast Fourier transform technique for solving Eq. ||| serves as a check on the alternate method 
represented in Eq. || which involves a set of coupled integral equations which can be solved 
by vector integration. Before presenting results (similar work can be found in refernece&^ii) 
for the solution of the gap equation, we present, for completeness, the expression for the 
penetration depth A^ 2 (T) which is a tensor. Results for this quantity will be given in the 
next section. The expression for A7?(T) is: 

~df(E k ) df(e k y 



\ -2 Cfe k 



(9) 



where the ^ are related to Fermi velocities and the overall constant in Eq. [|is of no interest 
here. 

In Fig. [I|, we show a phase diagram for the gap at zero temperature as a function of filling 
n (see Eq. ||) for a case t = lOOmeV, B = 0.45, T = and g XyV = 0.751 A small orthorhombic 
distortion 5 = 0.1 is also included in the calculations. Throughout the phase diagram, only 
the two components d x 2_ y 2 oc (cos^) — cos(k y )) and s x 2 +y 2 oc (cos^) + cos(k y )) appear in 
the gap. The amplitude multiplying these two functions are denoted by solid circles (d x 2_ y 2 
part) and solid square (s x 2 +y 2 part). The relative phase between the two basis functions 
is also given. We note that in the region down to approximately n = 0.45, the solution at 
zero temperature T = is mainly (i-wave with a small admixture of extended s-wave. The 
relative phase between these two components is zero throughout this region which means 
that the zero line nodes of the gap will be displaced somewhat off the main diagonals. This 
is seen in the angular resolved photoemission data of Ding et al.(13) although our model 
includes a single plane and nodes are found only on one side of the diagonals and not on 
both sides as seen in the experiments. This can be understood, however, quite naturally if 
a bilayer model were use (flEl to model Bi 2 Sr2CaCu 2 8 instead of the single plane model 
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used here. The small admixture of s x 2 +y 2 just described is due entirely to our use of a small 
orthorhombic distortion in the electronic dispersion0 and the amplitude of this component 
would be zero in a pure tetragonal case so that we would get instead, in this case, a region 
of pure d x 2_ y 2 in our phase diagram. 

Next, we note in Fig. [I] that below approximately n = 0.45 up to about 0.28, we have a 
region of admixture of extended s-wave and d x 2_ y 2 coexistence but with a finite relative phase 
of (p. For = 7r/2, the admixture is roughly of equal proportion while a small relative phase 
corresponds to a small amount of extended s and a phase near tt a small amount of d x 2_ y 2. 
Below n w 0.28, we recover a region of zero relative phase but with the two amplitudes 
having opposite signs and the s x 2 +y 2 part dominating. If an orthorhombic distortion was not 
included in our calculation the phase would always be 7r/2 but the s x 2 +y 2 component would 
only set in at some critical hole concentration n c \ and the d x 2_ y 2 component would drop 
to zero sharply at some other smaller critical concentration n C 2 with n c i and n C 2 dependent 
on the parameters of the theory such as details of the pairing potential and the underlying 
band structure. 

The region of finite phase (not zero or tt) in Fig. [I| is very interesting and was investigated 
further. For a fixed value of filling n taken to be 0.36, we show in Fig. Q our results for 
the two components of the g function of temperature. At high temperatures near 

the critical temperature T c , we have a sharply rising d x 2_ y 2 (solid circles) component as T 
is reduced and a small s x 2 +y 2 (solid squares) component. There is zero phase between these 
two components. But abruptly at T « 36K to be compared with a critical temperature 
of approximately 68K the phase switches on and the s x 2 +y 2 component starts rising quite 
rapidly. Our results are in qualitative agreement on these points with the angular resolved 
photoemission data of Ma et al.0 Of course our phase diagrams only roughly approximate 
the observed doping dependance. In view of this we make no attempt to fit the data since 
our model is quite simplified and, in particular, does not include the bilayer nature of the 
Cu02 planes in Bi 2 Sr 2 CaCu208, and so in the high temperature regime, our gap nodes 
fall only on one side of the diagonals. At lower temperature, once there is a finite relative 
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phase between the two components, the line node in the gap at the Fermi surface is lost, of 
course, because it is no longer real and both real and imaginary components would need to 
vanish simultaneously. For a d x 2_ y 2 and s x 2 +y 2 mixture, this can happen only at four points 
(±7r/2, ±7t/2) which, in general, will not fall on the Fermi surface. This can be seen in much 
more detail in Fig. |3|, but before doing so, we point out the insert in Fig. |2| in which we show 
the unnormalized penetration depth in the two principles in plane directions as a function 
of temperature. Since 8 ^ in our calculations, \ xx and X yy will be different although, 
as can be seen from the figure, their temperature dependence tracks each other fairly well 
and both show a clear change in slope at T = 36K, the point at which the relative phase <fr 
becomes finite. 

In Fig. 0, we show contour plots of the gap in the two dimensional CuC>2 plane Brillouin 
zone as well as a plot of the Fermi surface. Four frames are shown. Two correspond to 
the region above the second critical temperature T C 2 ~ 36K and two below. Above T C 2, 
the phase between the two components of the gap is zero and so the order parameter is 
real and its value can be plotted as a function of k x and k y . We see in Fig. |3] (frame b) 
positive and negative contours for A(T = 40K) and that the zero contour is not a straight 
line along the diagonal as it would be for pure d x 2_ y 2 wave. Instead, it is displaced below 
the diagonal around (ir, n) and above it around (0,0). Only at (tt/2,tt/2) is it exactly on 
the diagonal. Referring again to Fig. it is clear that the amount the node in the gap at 
the Fermi surface will have moved off the main diagonals will depend on temperature. In 
frame (d) of Fig. |3|, we show |A(T = 40K)| rather than its actual value. This graph differs 
from frame (b) only in as much as all contours are for positive values of the gap. It is to be 
used for easier comparison with the data of frame (a) which is for |A(T = 30K)|. In this 
case, the gap is complex and there is no zero contour for |A(T = 30K)|. The amount of 
change in the contours as compared to the previous frame (d) is not large, however, since 
the temperature has not been changed by much. The final set of contours shown in frame 
(c) are for low temperature and now these are quite different. We note, in particular, that 
|A(T = IK) | has a value near 6meV at the Fermi surface in the direction (n, it) and about 
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16meV for the (0, n) direction in rough agreement with the experimental situation found in 
Bi 2 Sr 2 CaCu 2 8 .(17) 

Using a simple tight binding model for the electron dispersion relation and a pairing 
potential involving a nearest neighbour attraction, we have calculated a phase diagram for 
the gap as a function of filling including in the calculation a small orthorhombic distortion. 
At high filling, the gap is found to be mainly of d x 2_ y 2 symmetry with a small admixture of 
extended s- wave (s x 2 +y 2) less than 10% in our work. The phase between the two components 
at zero temperature is zero so that in this region of phase space there are line nodes and 
the zeros of the gap in the Brillouin zone are slightly displaced off the diagonal. The 
degree of displacement increases with decreasing filling. In our work, the s x 2 +y 2 component 
is completely due to the orthorhombic distortion and disappears in a tetragonal system. 
At some critical filling n c i, a phase develops between d x 2_ y 2 and s x 2 +y 2 components. As 
the phase grows rapidly with decreasing filling, the amount of d x 2_ y 2 component decreases 
rapidly and the amount of s x 2 +y 2 increases. At some lower critical filling n c2 not so different 
from n c i, the phase becomes n and remains at this value for n < n c2 in which case the state is 
mainly s x 2 +y 2 with a smaller admixture of d x 2_ y 2. In a tetragonal system, the concentrations 
n c i and n c2 would be between the boundary of pure d x 2_ y 2 and mixture d + is and between 
d + is region and pure s wave respectively with the phase in the mixed region exactly ir/2. 

In the region of finite phase for the zero temperature gap, the temperature variation of 
the gap is complicated. At high temperature near T c (equal to 68K for the case considered), 
the gap is real and is an admixture of mainly d x 2_ y 2 with a small amount of s x 2 +y 2 so that the 
zero gap contours are displaced from the main diagonals. At some critical temperature equal 
to 36K in our example, a finite phase develops between the s and d components. While the 
d x 2_ y 2 component changes little at this temperature and shows an overall temperature de- 
pendence that is similar to that of an isotropic superconductor, the s x 2 +y 2 abruptly increases 
with decreasing temperature as observed in the angular resolved photoemission experiments 
of Ma et al.(17) In this temperature region, the gap shows no zeros on the Fermi surface 
and in our model ends up at T = with s x 2 +y 2 part roughly half the d x 2_ y 2 part. These 
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ratios depend on the parameters used of course. Calculation of the temperature dependent 
penetration depth show a clear change in slope of both the x and y components of the in 
plane penetration depth, \~£ and \~y, at T = 36K. 

Research supported in part by the Natural Sciences and Engineering Research Council of 
Canada (nserc) and by the Canadian Institute for Advanced Research (ciar). We thank 
R Soininen for discussions and interest. 
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FIGURES 

FIG. 1. Phase diagram for the gap at zero temperature as a function of hole concentration with 
n = 1 corresponding to half filling. The parameters used in the calculations are t = lOOmeV for the 
nearest neighbour hopping parameter, B = 0.45 for next nearest neighbour hopping, 5 = 0.1 for the 
orthorhombic distortion in the dispersion curves and g X;V = 0.75t for the nearest neighbour pairing 
potential. The solid circles give the amplitude of the d x i_ y 2 component and the solid squares the 
extended s x 2 +y 2 component. The gap in meV is given by the left hand side vertical axis while the 
phase between the two components is given by the open circles and the scale is shown on the right 
hand side vertical axis. It varies from to tt. 

FIG. 2. The temperature dependence (T) of the gap for a case with filling n = 0.36 which 
falls in the region of finite phase (j) and complex gap in the phase diagram of Fig. 1. The other 
parameters are as given in Fig. 1. The critical temperature is 68K. The amplitude of the d x 2_ y 2 
component is given by the solid circles and that of the s x 2 +y 2 by the solid squares. The phase 
which sets in abruptly at T = 36K is shown as the open circles. Again, the scale for the gap is on 
the left hand side in meV and for the phase on the right hand vertical axis in degrees from (0 to 
90). The inset shows the temperature dependance of the x and y penetration depths. 

FIG. 3. Contours of the gap as a function of momentum in the first Brillouin zone of the Cu02 
plane. Also shown is the Fermi line for filling n = 0.36. Because of the orthorhombic distortion 
included in our work through the parameter 8 in the dispersion curve, the Fermi surface does not 
have tetragonal symmetry and is open in the horizontal direction and crosses the gap contours. 
Frame (a) is for |A(T = 30K)|, frame (b) for A(T = 40K), frame (c) for |A(T = IK) | and frame (d) 
for |A(T = 40K)|. Note that the value of the basis functions at (0,7r) is two times the amplitude 
in Fig. 2. 
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